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A novel approach, the fermion-spin transformation to implement the charge-spin 
separation, is developed to study the low-dimensional t-J model. In this approach, 
the charge and spin degrees of freedom of the physical electron are separated, and 
the charge degree of freedom is represented by a spinless fermion while the spin 
degree of freedom is represented by a hard-core boson. The on-site local constraint 
for single occupancy is satisfied even in the mean-field approximation and the sum 
rule for the physical electron is obeyed. This approach can be applied to both one 
and two-dimensional systems. In the one-dimensional case, the spinon as well as the 
physical electron behaves like Luttinger liquids. We have obtained a gapless charge 
and spin excitation spectrum, a good ground state energy, and a reasonable electron- 
momentum distribution within the mean-field approximation. The correct exponents 
of the correlation functions and momentum distribution are also obtained if the 
squeezing effect and rearrangement of the spin configurations are taken into account. 
In the two-dimensional case, within the mean-field approximation the magnetized fiux 



state with a gap in the spinon spectrum has the lowest energy at half-filhng. The 
antiferromagnetic long-range order is destroyed by hole doping of the order ~ 10-=- 15 
% for t/J = 3-h 5 and a disordered flux state with gapless spinon spectrum becomes 
stable. The calculated specific heat is roughly consistent with observed results on 
copper oxide superconductors. The possible phase separation is also discussed at the 
mean-field level. 



PACS numbers: 71.45. -d, 75.10. Jm 
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I. INTRODUCTION 



The large U Hubbard model and its equivalent, the t — J model are prototypes to study 
the strong correlation effects in solids, especially in connection with the high supercon- 
ductivity |l]J3Q]. The central issue under debate is whether the non- Fermi liquid behavior, 
showing up as charge-spin separation and vanishing of the quasi-particle residue, inherent 
to the one-dimensional (ID) Hubbard model is also true for two-dimensional (2D) models, 
as conjectured by Anderson [Q. 

In ID, the exact Bethe-ansatz solutions [^,^ are available for the t-J model in the limit 
J/t and J/t = 2. The Hubbard model and the t-J model in the small J limit scale 
to the Luttinger model f^J^J^. Using Lieb-Wu's exact wave function, Ogata and Shiba |^ 
have shown the existence of an electron Fermi surface as well as a singular behavior at 
k ^ kp and ~ S/ci? in the electron- momentum distribution function. Moreover, Yokoyama 
and Ogata 0, and Assaad and Wiirtz have studied the ID t-J model using the exact 
diagonalization of small systems and quantum Monte Carlo methods, respectively, and their 
results show that the t-J model behaves like Luttinger liquids for low values of J/t, and 
undergoes phase separation at large values of J/t. Hellberg and Mele [|lT| came to the same 
conclusion by using the Jastrow variational wave function. Thus the typical behavior of the 
Luttinger liquid |T2| in ID, i.e., the absence of quasi-particle propagation and charge-spin 
separation, has been demonstrated explicitly for the t-J model in the small J limit. 



There are no exact solutions available in 2D. The variational calculations ||13[ seems to 
support Anderson's conjecture. The quantum Monte Carlo simulations gave some hint at 



vanishing of the quasi-particle residue in the thermodynamic limit ||14[. However, this result 
is not conclusive because of the " fermion minus sign" problem in the Monte Carlo technique 
and contrary results in exact diagonalization of clusters [|15i as well as analytic treatments 



of the single- hole problem [If: 



The crucial requirement |jT^ for the t-J model (and the large U Hubbard model) is to 
impose the single occupancy constraint ScrC'iLC'jcr < 1- An intuitively appealing approach 
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to implement this constraint and the charge-spin separation scheme is the slave particle 
formalism [0,0, where the electron operator is decomposed as Cjo-=a|/jo- with a\ as the 
slave boson and fi^ as fermion and the local constraint Y.a fLfi^r + '^I'^j = 1, or vice versa, 
i.e., a| as fermion and fia as boson. Due to the constraint, these particles are also coupled 
by a strong gauge field [|TU]|, allowed by this slave-particle representation. In the mean field 



approximation (MFA) the spin (spinous) and charge (holons) degrees of freedom are fully 
separated. However, there are a number of difficulties in this approach. First of all, in 
the slave boson version, the antiferromagnetic correlation is absent for zero doping, so the 
ground state energy in the 2D case is high compared with the numerical estimate for small 
clusters, and the Marshall sign rule ||2^ is not obeyed pT] , p2| . Alternatively, in the slave 



fermion approach, the ground state is antiferromagnetic for the undoped case and persists 



until very high doping (~ 60 %) ||23|. The large Fermi surface of spinous, present in the slave 
boson approach is absent there. Moreover, if we, following the common practice, let ficr keep 
track of the spin, while keep track of the charge, satisfying the sum rules: 6 =< a|aj > 
and 1 — 6 = J^a < fLfi(T where 6 is the hole doping concentration, we find that the 
sum rule for the physical electron J2a < Cj^Ci^ >= 1 — 5 is not satisfied for both versions. 
This expectation value is (1 — 5)^ in the slave fermion representation, and 1 — 5^ in the slave 
boson version. Since the total number of particles does not depend on the interactions, this 
difficulty will persist even beyond the MFA, so long as the spinon and holon expectation 



value decoupling is assumed. Furthermore, we have also shown p4[ that the overall electron 
distribution does not have the appropriate Fermi surface within this scheme even for the ID 
case. These are intrinsic difficulties of this decoupling scheme. 

In this paper we develop a new scheme, the fermion-spin transformation, to implement 
the local constraint and the charge-spin separation. In this scheme the charge degree of 
freedom is represented by a spinless fermion, while the spin degree of freedom is represented 
by a hard-core boson in terms of Pauli operators (with a projection operator to be specified 
later). Using this representation the local constraint is satisfied in the decoupling scheme in 
contrast with the existing slave particle approach |1I8|J3[], where the local constraint is replaced 



by a global one. As a consequence, the sum rule for the physical electron is obeyed. Moreover, 
the hard-core bosons can be expressed in terms of spinless fermions via the Jordan- Wigner 
transformation in ID and its generalization in 2D |^. This is an efficient calculation 
scheme which can provide very good results even at the mean-field level. 

Here we summarize our main results. In ID we can integrate out the spinless charge field 
(holons) and obtain an effective Hamiltonian for an interacting spinon field which behaves 
like a Luttinger liquid. Hence the physical electron, as a convolution of spinon and holon, 
also behaves like a Luttinger liquid in consistency with the exact solution Moreover, we 
obtain a gapless spectrum for both holons and spinous at the mean-field level which is not 
true in the slave fermion approach. The ground state energy at and away from half-filling is 
in good agreement with exact results. By going beyond the MFA to include the " squeezinc}'' 
effect and rearrangement of the spin configurations due to the hole presence, we obtain not 
only correct exponents of correlation functions and momentum distribution at the Fermi 
surface but also a reasonable global distribution function. In 2D we have considered various 
phases at and away from half-filling in the MFA. The magnetized flux state with a gap in 
the spinon spectrum has the lowest energy at half-filling. The antiferromagnetic long-range 
order (AFLRO) fades away by hole doping of the order lO-i-15% for t/J = 3^5in contrast to 
the Schwinger boson approach where the AFLRO is destroyed at 62% doping [^. Beyond 



the critical concentration, a disordered fiux phase with gapless spectrum becomes stable. 
We have also calculated the specific heat and considered the phase separation issue. The 
results are consistent with experiments and numerical simulations, respectively. 

The rest of the paper is organized as follows: In Sec. II, we explain in detail the fermion- 
spin transformation which is exact in the single occupancy Hilbert space, if a projection 
operator is introduced to remove the extra degrees of freedom. We also estimate the errors 
introduced by the MFA. In Sec. HI we apply the proposed scheme to ID t-J model within 
the MFA. The main results obtained have been mentioned above. In Sec. IV we calculate 
the correlation functions and momentum distribution by introducing two "string" operators 
which take care of the " squeezinc}^ effect and rearrangement of the spin configurations. The 



exponents thus obtained agree with the exact results. The apphcations to 2D at the mean- 
field level are described in Sec. V. Finally, in the concluding Section we make some further 
remarks to explain our current understanding why this simple transformation works so well 
and outline some open problems. 

II. FERMION-SPIN TRANSFORMATION TO IMPLEMENT THE 
CHARGE-SPIN SEPARATION 

A. Model, Constraints and Sum Rules 

We start from the t-J model which can be written as 

H=-tY: (ClCj. + h.c.) - /i E ClC,„ + JY.^.- S„ (1) 

{ij)a- icr (ij) 

where Cj^ (Cia) are the electron creation (annihilation) operators, Sj = CjaCi/2 spin op- 
erators with cr as Pauli matrices, fi the chemical potential. The summation (ij) is carried 
over nearest neighbour nonrepeated bonds. 

The Hamiltonian (1) is defined in a restricted Hilbert space without double electron 
occupancy. There are two ways to implement this requirement: either to solve (1) com- 
bined with a nonholonomic constraint J2a Cj^Cit^ < 1 or to introduce constrained fermion 
operators ||2^, replacing Cjo- by Ci„ = Cio-(l — rii^), where rii^ = Cj^Ci^. We will use both 
representations in this article. 

The constrained operators Ci„ satisfy the following relations 

Y: cia^ = y: cia^ii - n,,), (Y etc..) = 1-5, (2) 

where the latter equation is a sum rule for the electron at the hole doping concentration 
5, and (• ■ ■) means thermodynamical average. The on-site anticommutation relation of the 
constrained electron operator Ci^ is 

Y{c^.. cl} = 2 - E ctQ., (E{a., Cl}) = 1 + 5, (3) 

cr (7 (T 
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which gives rise to a sum rule for the spectral function Accr{k,w) 

j:r ^A,^{k,w)=i+6. (4) 

J-oo ZTT 

Of course, this value is less than 2 since 1 — 5 states are pushed to infinity as f/ — oo in 
deriving the t-J model. Hence the electron spectral function Accr{k,w) describes only the 
lower Hubbard band. Eqs. (2) and (4) are exact sum rules for the t-J model, and they must 
be preserved in adequate treatments. 

B. CPi Hard-Core Boson 

The decoupling of the charge and spin degrees of freedom for the physical electron is 
undoubtedly correct in the ID t-J model but the situation is still not clear in 2D. In this 
paper, we presume that the decoupling of charge and spin degrees of freedom for physical 
electron is also valid for the 2D t-J model, and propose a new scheme to decouple the charge 
and spin degrees of freedom. 

To motivate this transformation, we start from the no-double occupancy local constraint 
J2aCjcrCia < 1. Suppose Cia = h\hia^ wlth the spinless fermion operator hi keeping track 
of the charge (holon) while the operator keeping track of the spin (spinon), then this 
on-site local constraint can be rewritten as Z^o- ^i^l^L^io- ^ 1- Since the electron obeys the 
Fermi statistics, the operator hi^j must be a boson when the operator hi has been assigned a 
fermion character in this electron decoupling scheme. If bosons are subject to the condition 
Y.a^lj'itT = 1) the on-site electron local constraint 

= E ClC,, = hih\ = 1 - h% < 1 (5) 

cr 

is exactly satisfied, where rih = h\hi is the spinless holon number at site z, equal to 1 or 0. 
This way the nonholonomic on-site electron constraint is converted into a holonomic boson 
constraint. 

We should note that so long as h\hi = 1, I]<t C'io-C'io- = 0, no matter what is the value 
nb = Hcj^la^ia- However, the choice = 1 is convenient, because it also guarantees the 
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condition Uc = 1, when Uh = 0. This decouphng scheme, so called CP^ representation was 
proposed by Wiegmann |^ and was used in Ref. 29. The constraint ric = 1 means the 
presence of one boson (spin- up or -down) on each site, i.e., we assign a "spin" even to an 
empty site. This will not affect the physical expectation values, because the hole number 
expectation (ri^) will remove the spurious effects. Nevertheless, the extra degrees of freedom 
will affect the partition function and thermodynamical quantities. In the next subsection we 
will define a projection operator to cure this defect. As a result, the commutation relations 
and sum rules (2)- (4) will be satisfied exactly. 

Now we explore further the properties of the CP^ bosons. First note if we restrict the 
boson occupation number to or 1, the infinite-dimensional Fock space for bosons become 
two-dimensional, where we can choose the following representation 

/O 0\ /O 1\ 

which are nothing but spin-loweringS"" and spin-raising operators for = 1/2 and satisfy 
the hard-core constraints bb = Wb"^ = 0. 

Moreover, if we request that | and | hard-core bosons satisfy the CP^ condition = 1, 
the 2x2 representation space becomes two-dimensional. Assume , are singly 

occupied and empty spin-up, while , are singly occupied and empty spin-down 

VIA VoA 

states, respectively. Due to the constraint bj-^bi^ + = of 4 possible states as 

direct products, only two, namely , and are allowed. Thus we 

/O 0\ /O 1 

can ignore the spin label in the state and represent 6| as , and bi as 

VI oj \0 0, 

in the reduced two-dimensional representation space. Of course, all the hard-core boson 

conditions, i.e., biab\^ + b\„biu = 1, b\„b\^ = bi^bia = 0, (without summation over a), are 

satisfied. As a result, b^ and bi are identified with the spin lowering S~ and raising 

operators, respectively, while the boson occupation space is identified with the spin 1/2 

representation space. 
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To sum up, as solutions of the single occupancy constraint J2a C\„Cia- < 1 under CP^ 
convention J^ablabia = 1, we find the following fermion-spin transformation 

Cii = h\Sl , Ci[ = h\Sf (7) 

in terms of which the t-J Hamiltonian (1) can be rewritten as 

H = -tY. hih]{StSj + SiS+) + h.c. 
-^lY.hih\ + Jm^- h\h,)S, ■ S,(l - h]h^), (8) 

where Si is the pseudo-spin operator at site i which can be expressed as CP^ hard-core 
bosons and is different from the electron spin operator in Eq. (1). We would emphasize 
that the present CP^ hard-core boson representation of spin operators is different from the 
CP^ boson representation of spin operator used before p8|,p9|. In their approach, the spin 



degree of freedom is represented by ordinary boson operators, while, in the present scheme, 
the hard-core boson operator bier behaves as a fermion on the same site, and as a boson on 
different sites. 



C. Projection Operator 

In the local representation the restricted Hilbert space of no-double occupancy consists 
of three states, |0), | t), | |), while in the fermion-spin transformation presented in the 
previous Subsection there are four states \hole) ® \spin), namely |1,T)) ll;!)) |0,T)) and 
|0, i), where 1 or means hole occupation. We can introduce a projection operator P to 
remove the extra degrees of freedom. The matrix elements of this operator can be defined 
as 

P,^ = \K){al (9) 

where \k) is one of the bases of the physical states, while |a) is one of the bases in the space 
\hole) (g) \spin). Since the space dimensions of \k) and |a) are different, the usual relations 
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for the projection operator P"^ = P = are not satisfied. Using this operator, one can 
define the electron operators in the constrained space as 



Cjt — Pih\Si P} , Cii — Pih\Sf P} 



i] — PihiS^ Pj , (7/| = PihiS^ P/, 



(10) 



where Pi is the projection operator for the site P} is the hermitian conjugate of Pj. Making 

/O 1\ /O 0\ 

use of the matrix representation of the holon operator /ij = , /ij = , we can 

Vo 0/ \i oj 

write down exphcitly all these operators in matrix form (see Appendix). In particular, the 
constrained electron operators in the basis of the physical states |0), | t), | |) can be written 

as 

/O 1 0\ /O 0\ 

0, = 1 



a 



it 



il 



yo 0/ 

/o i\ 



VO 0/ 



yo 0/ 

/o o\ 



VI 0/ 



(11) 



which are nothing but the Hubbard X operators Xof, Xq^, etc. |30]. It is then straightforward 
to check that 



(12) 



(13) 



where 



(I 0\ 







is the hole number operator. Taking expectation values of Eqs. 



Vo 0/ 

(12) and (13), one sees immediately that the sum rules (2) and (4) are exactly satisfied. Thus 
we have shown that the fermion-spin transformation defined with an additional projection 
operator P satisfies exactly the no-double occupancy constraint and all sum rules,'i.e., they 
are an exact mapping. 
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However, the projection operator P is cumbersome to handle and in many cases, for 
example in the mean- field treatment of Sections III-V, we will drop it. Now let us see 
which of these properties are still preserved and what kind of errors we are committing in 
such approximate treatments. First of all, the local constraints are exactly obeyed even 
in the MFA. Secondly, those expectation values of electron operators, including spin-spin 
correlation functions, which should vanish, actually do not appear due to the presence of 
the holon number operator Uh = h\hi. Furthermore, as we will see later, the sum rules for 
the physical electron are also satisfied in MFA. By adding an extra spin degree of freedom 
to an empty site we are making errors of the order 5 in counting the number of states, 
which is negligible if 5 — 0. For comparison we should note that in the usual slave-particle 
approach |18|, the local constraint is explicitly replaced by a global constraint in MFA, 



and therefore the representation space is much larger than the representation space for the 



physical electron, which leads to some unphysical results [Il7| , p4| . From this point of view, 
our treatment of constraints for the physical electron is much better than the slave-particle 
approach, and therefore we believe the mean-field result based on the fermion-spin approach 
even without projection operator should be better than those obtained by using the slave- 
particle mean-field theory. This is indeed confirmed by the mean-field calculations presented 
in the following Sections. 

We note that a similar transformation has been discussed in Ref. 31, but these authors 
did not stick to the single-occupancy constraints in their actual calculations. 

III. THE MEAN FIELD THEORY IN ID WITHIN THE FERMION-SPIN 

APPROACH 



Since an exact solution ||3^ for the Hubbard model (hence for the J/t ^ limit of 
the t-J model) is available in ID, it is important to confront any approximate treatment 
with this solution. In this Section we consider the ID t-J model, using the fermion-spin 
transformation described in the previous Section at the mean-field level, neglecting the effects 
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of the projection operator P. As mentioned in Introduction, the ID t-J model exhibits a 
Luttinger hquid behavior, including charge-spin separation, vanishing of the quasi-particle- 
residue, etc. We should also mention that both spin and charge excitations are gapless in 
ID [0. In the standard slave-particle approach, many of these properties are not preserved. 
For example, in the slave-boson case there is a Bose condensation at the mean-field 
level which leads to a Fermi-liquid behavior. On the other hand, there is a gap in the spin 



excitation spectrum in ID within the slave- fermion framework, even at half-filling ||23|. Now 
consider the mean- field results in the fermion- spin approach. 

A. Luttinger liquid behavior 

In the fermion-spin representation, the ID t-J model may be written as 
H = -tY,aUi+i{StS~+i + + h.c. 

i 

-/i^ajoi + J^(aJai)(Si ■ Sj)(a]aj), (14) 

where, for the convenience of the following discussion, we have introduced the particle op- 
erators tti and aj defined as 



Qi = hi, a| = hi. (15) 

In the ID case, the hard-core bosons and can be mapped exactly onto spinless fermions 
using the Jordan- Wigner transformation 

St = fle'''^^<^^^^ (16) 

5'- = /.e-*"^'<'^'^', (17) 

Sf = fIf^-l, (18) 

where /« is the spinless fermion operator. Substituting Eqs. (16)-(18) into Eq. (14), the t-J 
Hamiltonian (14) can be expressed as 
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H = -tY^ial^i+i + al+i^difi fi+i + fi+ifi) - 

i i 

+jJ2ala,[^{flf) + f]f,) + (///, - \)U]f, - \)]a]a,. (19) 

We can now employ the path-integral representation in which the Lagrangian L and the 
partition function Z of the t-J model in the imaginary time r can be expressed as 

L = Y, a\drai + fldrh + H, (20) 

i i 

Z = J DaDa^DfDfe-f'^^^'-^\ (21) 

Integrating out the spinless charge field of the t-J model, we obtain an interacting spinon 
system, which like any interacting fermion systems in ID, is described by a Luttinger liquid 
theory |T2| , p^ . Moreover, the physical electron as a convolution of spinon and holon also 



behaves like a Luttinger liquid, which means that the electron wave function renormalization 
constant Z = in the ID t-J model. 

In the path-integral representation, one can introduce a 5'f/(2)-invariant Hubbard- 
Stratonovich transformation to decouple the Lagrangian (20) by using the following auxiliary 
fields 

Xi,i+V ~ ^t^i+r)^ (22) 
4'i,i+ri = O^^O^i+rii (23) 

where r] = ±1. Note the auxiliary field Xi,i+ri is a boson type field. MFA to the t-J model 
(19) is just the saddle point solution of the Lagrangian (20), i.e., the auxiliary fields Xi,i+v 
and (pi^i+n are replaced by their mean-field values Xi^+v = X and 0j,i+r, = 4>, respectively, and 
the Hamiltonian (19) can be diagonalized as 

H = Y.e{k)alak + T.^{^)flh + 4iVtx0 + 2iVJ[(l - 5f - ^'■\x\ (24) 

k k 

where N is the number of sites, and 

13 



e{k) = —4:txcos{k) — fi, 



(25) 



coik) = [((1 - 5f - 02)(1 - 2x) - j<P]{2J)cosik), (26) 

while the self-consistent equations for the order parameters x ^inci can be obtained by 
minimizing the free energy. We can now proceed to a brief discussion of the results in MFA. 



B. Groud-State Properties at Half-Filling 

At half-filling, there are no charge degrees of freedom, and the t-J model (24) reduces to 
the antiferromagnetic Heisenberg model in the fermion representation 

Hj = Y.Mk)flfk + 2NJx'. (27) 

k 

In this case the order parameter x can be evaluated to be x = — l/^r, and we obtain a gapless 
spinon spectrum, 

a;o(A;) = (l + -)(2J)cos(A;), (28) 

TT 

which is rather close to the exact result of the Heisenberg chain obtained by using the Bethe- 



ansatz method |^ uJo{k) = {n /2){2J)cos{k). Correspondingly, the spinon spectrum near the 
spinon Fermi surface [k = ±7r/2) is linear with velocity = (1 + 2/tt){2J) = 1.6366(2J), 
which is also very close to the exact result of the Heisenberg chain obtained by Haldane 
pGl Vs = {ji /2){2J) = 1.5708(2J). The ground state energy of the Heisenberg model at 
temperature T = is Eq = — 0.4196(2J) which is only 5.3 percent higher than the exact 



Bethe-ansatz value ||] of = -(1/4 - ln2)(2J) = -0.4431(2J). 

Thus in the half-filled case, the spinon has a gapless spectrum, and the spinon ground 
state energy can be described adequately by the Jordan- Wigner transformation within MFA. 
This case has already been considered earlier [B^,P^. 
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C. Ground-State Energy away from Half-Filling 



Away from half-filling, a gapless liolon spectrum is obtained 



2 

e{k) = -2t-cos{k), (29) 
vr 

from which the holon Fermi velocity is given by 

2 2 
Vh = 2t-sin[(l - 6)tt] = 2t-sin(57r), (30) 

and the holon ground state energy at temperature T = is obtained as 

2f 2 

E, = -{-)i-)sm{5n). (31) 

TT 71 

All these quantities differ from the corresponding exact values |^ of the ID large U limit 
Hubbard model e{k) = —2tcos{k), Vh = 2tein(57r), and Eh = — (^)sin(57r), only by a factor 
2/7r. 

At the same time, the gapless spinon spectrum at finite dopings becomes 

uik) = [((1 - 6)^ - ^^^^^)(1 + -) + ^-sin(57r)](2J)cos(A;), (32) 



TT^ TT' 
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with the spinon ground state energy 



which is again very close to the exact result of the spinon ground state energy |^ of the 
ID large U limit Hubbard model over the entire doping range. Here Eq is the ground state 
energy at half-filling. Therefore the total ground state energy of the ID t-J model in the 
fermion-spin approach within MFA can be expressed as 

Eg = Eh + Eg 

= -(^)(2)sin(H - (1 - 5r[l - 4;fM^]0.4196(2J). (34) 

The ground state energy and thermodynamical quantities for the ID Hubbard model in the 
atomic limit have been calculated long time back |^|. Our result for the ground state energy 
is rather close to theirs. 
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D. Momentum Distribution 



Now we turn to discuss the global features of the electron momentum distribution within 
MFA. This distribution for physical electron is defined as 

n{k) = T.{ClC,^) = ^ E e^'^'^^-^-HclG^.). (35) 

<^ Inicr 

In the present fermion-spin approach, neglecting the effects due to the projection operator 
defined in Sec. II C, this distribution function can be rewritten as 

Im 

Using the Jordan- Wigner transformation (16)-(18), n{k) can be further expressed as 

Im 

+/ie"'"^^'<'^^'^'~^^<™^^^^VrL))- (37) 
Since e^'^^^if' = 1 - 2///^, in MFA, we obtain 



oo 



n(k) = 1-6+ y (-l)™+i(-)"+i^^^^^!^cos(mA;), (38) 

which is plotted (solid line) in figure 1 for doping 6 = 0.5. For comparison, the corresponding 
curves for the slave-boson (dashed line) and the slave-fermion (dot-dashed line) cases are also 
given. It is obvious that > 1 — 5 in the slave-boson case and < 1 — 5 in the slave-fermion 
case, which is very far from a should-be electron momentum distribution. The integrated 
area under the curve is equal to (1 — 6^) in the slave-boson case and is (1 — 6)"^ in the slave- 
fermion case, while the correct value should be 1 — 5 l^^. The solid curve corresponding to 



our transformation, is closer to the exact result The integrated area is correct and the 
shape looks like a reasonable momentum distribution, i.e., in some part the distribution is 
greater, while in other part it is less than 1 — 5. To get a more accurate result (including 
the correct location of the Fermi surface and a correct slope at it) one should go beyond the 
MFA to include the spinon-holon interactions as discussed in the next Section. 
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Thus away from half-filling, the spinons and holons are decoupled completely, with the 
holon behaving like a spinless fermion, while the spinon has the Jordan- Wigner form in 
ID. The gapless spectra for both holons and spinons, as well as the ground state energy can 
be described adequately within the fermion-spin approach even in MFA. 

IV. THE ASYMPTOTIC BEHAVIOR OF CORRELATION FUNCTIONS IN THE 

ID CASE 

A. Motivation 



Interacting ID electron systems generally behave like Luttinger liquids |12] where the 
electron correlation functions show a power-law decay with unusual exponents. These sys- 
tems exhibit an electron Fermi surface with a correct Luttinger volume but the momentum 
distribution function is singular at the Fermi surface, also with unusual exponents p[. These 



exponents depend on the interaction strength. Haldane ||12[ has shown that the character- 
istics of Luttinger liquids can be calculated using the bosonization techniques. To get more 
insight into the problem let us consider the Bethe ansatz wave function for the Hubbard 
model in the large U limit (hence for the t-J model in the small J limit), derived by Ogata 
and Shiba p 

^(a;i, ■ ■ -, xn) = det[exp{ikiXQj]^{yi, ■ ■ -, i/m), (39) 

where the determinant depends only on the coordinates xq. of particles, but not on their 
spins, while $(?/i, ■ ■ -,1/^/) is the Bethe ansatz wave function for an AF Heisenberg chain 



with yi, ■ ■ ■,yM as coordinates of down spins [35|. This asymptotic form can be interpreted 



as a complete separation of charge and spin degrees of freedom in some sense. In fact, the 
determinant describes spinless fermions (holons), whereas $ is the "spinon" wave function. 
Our fermion-spin transformation is, to some extent, an approximate second quantized version 
of this solution, with holon being represented by a spinless fermion and spin represented by 
a hard-core boson. 
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However, there is an important "detail" in the wave function (39), namely, the spin wave 
function $ is for a "squeezed^ Heisenberg chain, i.e., all holes should be removed. This will 
lead to rearrangement of spin configurations and far-reaching nonlocal effects. Therefore, 
spinous and holons are not completely independent, but interacting with each other rather 
strongly. As shown in Ref. 5, the correct exponents of correlation functions and an adequate 
description of the momentum distribution function (Fermi surface at kp, rather than 2kF 
with appropriate singular behavior) can be obtained only if these interaction effects are 



properly taken into account. Weng et al. have shown that the effects due to squeezing 
and rearrangement of spin configurations can be included by introducing a nonlocal " string" 
field. After doing that correct results for the correlation function exponents, etc., can be 
obtained for a ID t-J chain. In this Section we calculate these exponents within our fermion- 
spin approach, following their technique with some modifications. The nonlocal effects due 
to spinon-holon and spinon-spinon interactions will be included by introducing two string 
fields. After squeezing, there are no holes in the chain, so the additional degrees of freedom 
due to assigning "spin" to a hole site disappear, hence the projection operator introduced 
in Sec. II C is not needed for our purpose. 



B. String Operators 

Let us consider the the t-term in the t-J Hamiltonian within the fermion-spin repre- 
sentation. Following Ref. 29, the largest ho Ion kinetic energy may be obtained if the spin 
configurations are squeezed as 

Pi{alS^ P- Pi+iai+iSi_^_i + a-S'j P- Pi+iai+iS^_^i)P-_^_i = a-aj+i (40) 

for all sites i where a holon is present. However, after such squeezing, the spin configurations 
are not optimal to favor the spinon energy. Thus at the same time spin configurations must 
be rearranged into optimum configurations to provide the lowest spinon energy. These 
optimum spin configurations can be obtained by reversing the original spin polarization 
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for all sites on the left-hand side of site i. These processes of first squeezing the t — J 
chain and then rearranging the spin configurations, are shown schematically in figure 2. 
These nonlocal processes cannot be described by any formal perturbation theory, but can 
be taken into account approximately by introducing the string fields In our case, 

we introduce two string fields to describe the above physical processes, so the constrained 
electron operator Cia can be expressed within the fermion-spin approach as 

Ci^ = (e*-(^-E,>,'^J«Oa.)(g--E,<,'^J"^^-), (41) 

C,^ = (e*-(^+E,>,'^^Oa^)(g-E,<,'^^^^+), (42) 

which means that the spinless fermion and the hard-core boson Sf are replaced by cor- 
responding string operators as 

ai^a^e^'^^^^^'M'^'^'^ (43) 

^ S'fe^'"^'<«"^"S (44) 

where e^*'^^'<'"'"' is the string field for the spinon due to the presence of holons, describing 
the effects of rearrangements of spin configurations from -oo to site i when an electron is 
removed or added at site i. For the t — J chain, holons on the right-hand side of site i can 
feel some indirect effects of holons on the left-hand side of site i due to the rearrangements 
of spin configurations from -oo to site i when one electron was removed or added at site i. 
These indirect effects can be described by the string field e*'^*-^^^'>» for the holon. One 
can check easily that the anticommutation relations for the physical electron are preserved 
exactly in our case. 

C. Energy Spectrum for a Squeezed Chain 

After squeezing and rearranging the spin configurations, the t-J model can be written as 
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H = -t ^(a|a,+i + h.c.) -^iJ2^U^ + ^[(1 - - E ^^ ■ (45) 

where we have approximated the probabihty of the spin exchange process of the Heisenberg 
term as ajojajaj ^ {a\ai) {a^^aj) — {a\aj){a}jai) = (1 — 5)^ — 0^, and the lattice constant 
of the squeezed spin chain has become — 5), where a is the original lattice constant. 
Therefore the Fermi points of spinons are shifted from kp = ±7r/2 of the half- filled case to 
kp = ±7r(l — 6)/2 for the doped case. 

For the squeezed chain, the gapless holon spectrum, the holon Fermi velocity, and the 
holon ground state energy are 

e{k) = -2tcos{k), (46) 



Vh = 2tsin(57r), (47) 
Et = sin((57r), (48) 

71 

respectively, which are in full agreement with the corresponding exact values of the ID large 
U Hubbard model [^]. At the same time, the gapless spinon spectrum, spinon Fermi velocity, 
and the spinon ground state energy are given by 

uj{k) = [(1 - 6f - ^l^^M](2J)cos(A;), (49) 
V, = [(1 - - sin2(57r)/7r2](2J)cos(57r/2), (50) 



^.-[(l-^)^-^]i^o, (51) 



respectively, which are very close to the corresponding exact values of the ID large U limit 
Hubbard model Here Eq is the ground state energy at half-filling. Therefore, the total 
energy of the t-J model Eg = Et + Ej agrees quantitatively with the exact value of the 
ID large U limit Hubbard model in the entire doping range. It fully agrees with the ID 
results for the Hubbard model obtained earlier in the atomic limit [^. The gapless spinon 
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and holon spectra, and the ground state energy are all better than the mean-field results 
obtained in Section III, which indicates that the string fields take into account the spinon- 
holon interactions in the t-J model, and renormalize considerably the results obtained in 
MFA without string fields. 



D. Correlation Functions 

The spin-spin correlation function is defined as 

= {s^{t)sfm = {sj{t)s]m 

= {Sf{t)Sfm = \{St{t)Sr{0) + Sr{t)Sim. (52) 

In the doped case, we need to replace the operator Sf in Eq. (52) by Eq. (44) to account 
for the presence of holons. Thus substituting Eq. (44) into Eq. (52), we obtain the spin-spin 
correlation function of the t-J model as 

S{xi-Xj,t) = -[(S+(t)5,-(0))(e*"S^<-'^^(*)'^'(*)e-'"^'<^"'(°)"'(°)) 



4 

+ 



Following Haldane |12,29|, we apply the bosonization method to the free holon field, and 



obtain the following asymptotic form 



) r, (54) 

[{xi -XjY - {vhty]i 

where Vh is the holon Fermi velocity. Luther and Peschel p3 have mapped the ID Heisen- 



berg model into an interacting spinless fermion system by using the Jordan- Wigner p5 
transformation. Their work involves generalization of the Jordan- Wigner transformation to 
provide a representation for continuum spin operators. The asymptotic form of the spin- 
spin correlation function of the Heisenberg model can be then obtained by considering the 
spinon-spinon interactions. Following their calculation, we get 

(SmSiiO) + Srit)Sim ~ ^os[2fcHx.-x )] ^^^^ 
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where Vg is the spinon Fermi velocity. Combining Eq. (54) with Eq. (55), we obtain the 
asymptotic form of the spin-spin correlation function of the ID t-J model as 

Qf ,N 1 cos[2kF{xi - Xj)] 

S[Xi — Xj,t)r^ Y- — (56) 

which shows a power-law decay with exponent 3/2, in full agreememt with the numerical 
result of the ID large U limit Hubbard model obtained by Ogata and Shiba ||^. 

E. Momentum Distribution Function 

The electron momentum distribution function nik) is defined in Eq. (36) within the 
fermion-spin approach. To consider the squeezing effect and rearrangements of the spin 
configurations, we need to substitute Eqs. (43)- (44) into Eq. (36). Then n(/c) can be 
rewritten as 

where the factor Gi = e^^'^^~^i>i"'i"''-\ In what follows, we will drop this factor as in the 



previous calculations p9| , |24| , since it will only contribute a next-to- leading additional power- 
law decay in the asymptotic single electron Green's function, so one may neglect it if only 
interested in the leading contribution. 

The calculation for the global features of n{k) is similar to the case without string fields 
within MFA and the result is 

n{k) = l-5 + Ai{k) + A2{k) + A^{k) + ■ ■ -, 

2 

Ai{k) = (-)^sin(57r)cos(/i;), 

71 

9 1 A 

A^ik) = (-)3[-sin2(57r) - -sin(2(57r)]cos(2fc), 

TT TT 2 
1 1 OA 

A:i{k) = (-)^[-sin(57r)sin(257r) + -(1 )sin(357r)]cos(3A;). (58) 

TT TT 3 3 

The curve n{k) is plotted in figure 3 (solid line) for 6 = 0.25 in comparison with the result 
without string fields (dot-dashed line). We find some substantial improvement by including 
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the string fields. In particular, the location of the Fermi points was wrong in MFA without 
string fields, while it is correct {kp = |(1 — 5)) if they are included. In the same figure we 
have also plotted the n{k) curve obtained in the early treatment (dashed line) ||2^ for CP^ 
electron representation (without accounting for hard-core natuie of bosons). Obviously, the 
global features of this curve are not correct ( n{k) goes up again as k further increases). 
Nevertheless, the asymptotic form of the momentum distribution near kp, obtained in both 
approaches, is correct, namely, 

n{k) ~ const. — C\k — kp\»sgn{k — kp), (59) 

again in agreement with the exact numerical result . This singularity is washed out in the 
numerical calculations and does not show up in figure 3. 

To sum up we have confirmed that by introducing the string fields the spinon-holon 
interactions can be included to some extent which allows us to obtain correct exponents for 



the correlation functions and momentum distribution, as proposed earlier I^Q]. Moreover, 
the global features of the energy spectrum and the momentum distribution, found in our 
fermion-spin approach are correct in contrast with the previous approach |^ which could 
not provide such an adequate description. 

V. THE MEAN FIELD THEORY IN 2D WITHIN THE FERMION-SPIN 

APPROACH 

In this Section, we consider the 2D t-J model. Very soon after the discovery of oxide 
superconductors Anderson revived his idea of the resonating valence bond (RVB) to 
describe the short-range AF fluctuations in the 2D Hubbard model. Baskaran, Zou and 
Anderson developed a mean field theory of the RVB state. Later a number of other 
more elaborated mean field solutions have been discussed both at half filling and away from 
it, such as the flux phase the spiral phase [0,0 and the commensurate flux phase |^ 



which breaks the time-reversal symmetry and parity. The latter is related to the proposed 
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fractional statistics and anyon superconductivity For the early suggested flux phase 

P0| , the lowest energy solution is a state with | flux quantum (or phase vr) per plaquette 
which does not break the time-reversal symmetry and parity. Sheng, Su and Yu found that 



the TT-flux phase can coexist with s + id wave RVB state for small dopings MM . Lee and Feng 



2^ , and later Chen, Su, and Yu |^5| found that the magnetized RVB state with coexisting 



AF order and a d-wave RVB order parameter, has a gain in energy. Furthermore, Hsu |^ 
has shown that the magnetized flux state which is the coexistence of AF order with a flux 
state has a similar gain in energy. Recently, Wang and Ubbens and Lee [^, obtained 
the same result in a different framework. In this Section, we discuss the 2D t-J model along 
this line. 



A. Generalized Jordan- Wigner Transformation 



In the fermion-spin approach, the success of theory depends strongly on whether one can 
map the hard-core boson onto an appropriate fermion or boson representation. In the 2D 
case, the Jordan- Wigner transformation of the spin operators has been generalized by several 



authors [Q. In particular, Wang discussed the Heisenberg model in the MFA using 
this generalized Jordan- Wigner transformation. Using our fermion-spin transformation (7), 
this mean-fleld calculation can be easily generalized to the t-J model. The generalized 



Jordan- Wigner transformation may be written as p6|,|37 



(60) 



Si = fie 



sf = nu 



-idi 



1 

2' 



(61) 



(62) 



where is a spinless fermion, 6i = fi fiBa. In order to preserve the hard-core prop- 
erties of the spin operators, Bu should be such that e*'^" = —e"^^''. One possible way of 
materializing this equality is to set Bu = lmln{Zi — Zi), with a complex representation of 



24 



the lattice sites Zi = Xi + iYi. If the effects connected with the projection operator P are 
neglected, then the 2D t-J model can be mapped onto the fermion representation as 

H=-tJ2 h^h] iflfje'^'^-'^^ + h.c.) - /X ^ hM 

(*i> » 

+JY.h^h^^h,h][\{f}f,d^'^-'^^ + h.c.) + - !)(/]/, - \)]. (63) 

It has been shown p6| , |37| that the phase factor in the Hamiltonian (63) creates a gauge field, 
with the vector potential given by 

where ni = f} fi. On average, there is a vr-fiux attached to each spinless fermion /j, which 
is nothing but the Chern-Simons gauge field converting a hard-core boson into a spinless 
fermion. As mentioned earlier, the vr-fiux does not break the time-reversal symmetry and 
parity. 

The Hamiltonian (63) is obviously very complicated, so a more complete discussion about 
it beyond the MFA will be given elsewhere. In this Section, we only discuss some mean-field 
properties of the 2D t-J model within the fermion-spin approach. In the MFA, following 
the similar discussion of Laughlin et al. Mele and Wang the phase factor 



is absorbed by redefining a bond- dependent exchange parameter J^-. At half filling, there 
are many possible phases and the state with the lowest energy turned out to be the 
magnetized flux state with coexisting Neel order parameter {Sf) = (— 1)*M and orbital 
current order parameter {S^S^,,) = x, where z/ = ±x, ±y. The ground state energy per 
bond is Eq = — 0.33J, while the staggered magnetization is M = 0.389, which are rather 



close to the best numerical estimates Eq = — 0.3346J and M = 0.31, respectively. This 
state is completely equivalent to what was first discussed by Hsu who obtained the 
ground state energy Eq = — 0.331J and staggered magnetization M ^ 0.3 within a different 
theoretical framework. All these results show the accuracy of the mean-field result within 



the generalized Jordan- Wigner transformation [26|. 
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B. Mean-Field Theory away from Half-Filling 

Away from half-filling we need to introduce an additional holon particle-hole order pa- 
rameter = and the t-J model can be decoupled in the MFA as 

iJl) 3 W 

-2J,ffM E(-1)71/.- + 4A^^e//x' - SNtx<l> + ^NJ^ffM', (65) 

m 

where Jg// = J[(l — 5)^ — 0^], and N is number of lattice sites. In the magnetized flux 
phase, the t-J model can be diagonalized by using the Bogoliubov transformation to give 

^ = E i^kC^W + stPiPk) + E (EkAlA, - E.bIb,) 

k(red) k{red) 

+mJeffX^ - 9>NtX(t> + ^NJeffM^, (66) 

where (red) means the summation is carried only over the reduced Brillouin zone. The new 

operators a^, Pk, A^, and are related to h^, /if, f^, and fj^ by 

ht = ^{ak + Pk), h^ = ^ia,-P,), (67) 



and 



fk = '^'"t" {ukAk - VkBk), = UkBk + VkAk, (68) 

Jcos^ky + sin 



where 



2 1,^ 8JM 2 1- 8JM, 



t.^ = -(l + -^), vi^-{l-—^), (69) 



and the spin excitation spectrum is 



Ek = \l{^Jef}Mf + 4(Je// + 2t(t) - 2Jeffxy{cos^ky + siu^k,), (70) 
while the charge excitation spectrum is 

^fe = ±2tx7fc - 7fc = 2(cosA;^ cos/cj^). (71) 
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In obtaining the above results, we have considered two sublattices A and B with i & A and 
i + u E B. The free energy of the 2D t-J model can be obtained as 

F = -l ^ ln[(l + e-^^0(l + e-^^^)]-^ E ln[(l + e-^^'=)(l + e^^'')] 

^ k{red) ^ k(red) 

+ANJ,ffX^ -mtx<j) + UeffM\ (72) 

from which we find the self-consistent equations for the order parameters 4> by mini- 

mizing the free energy with respect to these parameters. 



C. Doping Dependence of the Staggered Magnetization 

At half-filling, the 2D t-J model reduces to a 2D AF Heisenberg Hamiltonian, and the 
result is the same as discussed by Hsu |^ and Wang [0. The spinon spectrum in the 
magnetized flux phase with AF long-range order is expressed as 



Ek = \/{SJMy + 4J2(i _ 2xy{cos'^ky + siii^k^), (73) 

where a gap appears due to the presence of AF staggered magnetization. We note that 
a gap in the spinon spectrum of the flux phase at half filling was suggested by Laughlin 
earlier However, this gap coming from the staggered magnetization M decreases very 



rapidly upon doping. In the MFA, M vanishes around doping 6 ~ 0.1 0.15 for t/ J = 3 -h 5 
which is plotted in figure 4(a). This result is in reasonably good agreement with experiments 
49| , |50[1 and Monte-Carlo simulations p2|. For comparison we note that the magnetization 



vanishes only at 6 ^ 0.62 in the MFA for the slave-fermion approach ||2^. At finite dopings 



but still within the magnetized fiux phase, we find a competition between the Neel order 
parameter M and the orbital current order parameter x, with M decreasing very rapidly 
(see figure 4(a)), and — x increasing very fast (see figure 4(b)) upon doping. In the small 
doping range, the holon particle-hole order parameter increases roughly linearly upon 
doping, and is almost independent of t/J, which is shown in figure 4(c). The t-J model 
is characterized by a competition between the kinetic energy (t) and the magnetic energy 
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(J). The magnetic energy J favors an AFLRO for the spins, whereas the kinetic energy 
t favors delocahzation of the holes and tends to destroy the spin AFLRO. Thus the rapid 
suppression of the AFLRO upon doping means that the present mean-field kinetic energy 
is better than those obtained in the slave-fermion approach and is closer to the optimal 
kinetic energy of the system. Beyond doping 5 ~ 0.1 -^ 0.15, corresponding to the range 
of the actual high-temperature superconductors, there is no AFLRO, but short range AF 
correlations persist and the spinous are in a disordered flux phase with a gapless spectrum. 



Ek = 2(Je// + 2t(f) - 2Je//x) ycos^^y + sin'K, (74) 
as conjectured by Anderson |T|. This spinon spectrum is similar to that of the flux state 



discussed also at the mean-field level by Affleck and Marston, and Kotliar pO |. 

The mean-field phase boundary between the magnetized flux and disordered flux states 
is, of course, at somewhat higher doping 6 than the value given by experiments and numer- 
ical simulations. In fact, the frustrations of spins can shift the mean-fleld phase boundary 



towards smaller doping S fSTj]. Thus we believe that the result will be even closer to experi- 



ments and numerical simulations by going beyond the MFA. 

D. Specific Heat 

The speciflc heat measurements on oxide superconductors have been made for many 
compounds by different researchers The descrepancics are mostly due to the difficulty 



in preparing and characterizing samples of the oxide compounds. Although the speciflc heat 
data for the superconducting compounds show considerable variations for samples measured 
by different groups, some qualitative features seem common to all the measurements. Hence 
a quantitative comparison between theory and experiment is still early, but the qualitative 
tendency of the speciflc heat in an adequate theoreticl description should be consistent with 
experiments. 

In the half fllled case, there are no charge degrees of freedom and the spinon speciflc 
heat has been considered by Wang Away from half-fllling, we are interested in the 
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doping range S = 0.1 ~ 0.2 where the superconductivity appears. In this doping range, we 
have shown that the magnetization M vanishes and the system is in a disordered flux state, 
where the internal energy of the system in the MFA is given by 



The numerical result is shown in flgure 5 at doping 6 = 0.2 for t/J=3 (solid line), t/J=5 
(dashed line), and the shape is similar to the experimental results [^^ . For T > 0.0005 J, 
the speciflc heat is found to increase with temperature T, which is also consistent with 
the experiments. Therefore our simple mean-fleld calculation provides a correct qualitative 
description of the speciflc heat for oxide compounds. 



The possible phase separation in the t-J model was proposed by Emery, Kivelson, and 
Lin [|^. They argued that the dilute holes in an antiferromagnet are unstable against phase 
separation into hole-rich and no-hole phases, and the transition from an ordered state to 
doped state is of flrst order. Later investigation |^ using high-temperature expansions 
shows that a line of the phase separation extends from J/t = 3.8 at zero fllling to J/t = 1.2 
near half-flUing, but for the range of parameter interesting to the copper oxide planes there is 
no evidence for phase separation. Within the mean-fleld theory of our fermion-spin approach, 
we flnd that the phase separation is robust for the t-J model, and the phase separation 
manifests itself at the mean-fleld level by a negative compressibility (slope of the chemical 
potential). The total energy Etotai and the chemical potential /i as a function of the doping 
6 for t/ J = 5 is plotted in flgure 6, which shows that the phase separation occurs roughly at 
dopings S < 6c = 0.08. In this doping range (5 < 6c), the compressibility of the t-J model 




(75) 



from which the speciflc heat can be obtained as 




(76) 



E. Phase Separation 
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is negative, and therefore the magnetized flux state with long-range order is unstable. The 
range of the phase separation will be reduced by considering the frustrations of spins . 



VI. SUMMARY AND DISCUSSIONS 

In this paper, we have developed a theoretical framework, fermion-spin approach, to 
study the low- dimensional t-J model. In this approach, the physical electron is decoupled 
as a product of a spinless fermion (holon) and hard-core boson (spinon). The main advantage 
of this approach is that the on-site local constraints of the t-J model or the large U limit 
Hubbard model can be treated exactly in analytical calculations. In this framework, we 
have shown that the holon behaves like a spinless fermion, while the spinon is neither boson 
nor fermion, but a hard-core boson, and the sum rule for the physical electron is obeyed. 
This is not true for the conventional slave-particle theories, where the spinon behaves like 
boson (slave-boson approach), or fermion (slave-fermion approach), and the sum rule for the 
physical electron is not obeyed within the decoupling scheme |^ . 



We have applied this approach to study the low-dimensional t-J model. 

In the ID case, we have obtained gapless spinon and holon spectra, a good ground 
state energy, and a reasonable electron-momentum distribution within the MFA. Thus the 
ground-state in the fermion-spin formalism is in some sense closer to the Bethe-ansatz Lieb- 
Wu's exact wave function of the ID large U limit Hubbard model than the slave-particle 
approach. It is shown that the spinon and consequently also the physical electron behave 
like Luttinger liquids. We have also obtained the correct asymptotic form of the spin-spin 
correlation functions as well as the electron-momentum distribution function of the ID t-J 
model within the fermion-spin approach by considering the string effects, with results in 
agreement with the exact numerical simulations of the ID large U limit Hubbard model 
obtained by Ogata and Shiba 0. 

The ID problem is a good testing ground where the charge and spin are truly separated 
(not in the literal sense of product of spinon and holon, but rather as independent collective 
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excitations) and the Fermi liquid theory fails to describe its physical properties. To our 
knowledge, neither the standard perturbation theory, nor the conventional slave particle 
approach is capable of handling this aspect. Our results of the MFA seem to hint that the 
fermion-spin approach has some potential to further explore. 

The results for 2D are also very encouraging. The magnetized flux phase with a gap 
proportional to the staggered magnetization, in the spinon excitation spectrum, has the 
lowest energy at half-filling in the MFA. However, this AFLRO fades away at 10 15% hole 
doping for t/ J = 3 5, beyond which a gapless flux state becomes stable. This result agrees 
with experiments and numerical simulations ||49| , |50|j2^ . It is essential that this mean- field 



result was obtained without any adjustable parameters. This means that the formalism 
itself is powerful enough to handle the frustration (delocalization) effect of the t-term in 
destroying the AFLRO favoured by the J-term. However, we should mention that the spin 
excitation spectrum at half-filling is not gapless in this approach within the MFA as it should 
be (spin waves). Both Laughlin's approach and Hus's treatment |^6|] suffer from the 



same weakness. Probably, it can be cured by including vertex corrections. 

As mentioned in Sec. II, this formalism doubly counts the empty site by assigning a 
"spin" to it. As shown there, this defect can be cured by introducing a projection operator 
to remove the extra degrees of freedom. However, in our mean-field treatments we have not 
taken this projection operator explicitly into account. The fact, that we still obtain very 
good results as summarized above, indicates that we are not making substantial errors by 
allowing this extra degree of freedom, at least in problems considered so far. 

A natural question is: What is the reason why this simple-minded transformation is 
so useful. To our present understanding, there are, at least, three reasons: (1) The local 
constraint is exactly satisfied even at the mean-field level. (2) The hard-core nature is kept 



in the calculation via the Jordan- Wigner transformation in ID [|25[ or its generalization in 
2D [^. (3) The representation of the hard-core boson in terms of spin raising and lowering 
operators is essential, because whenever a hole hops it gives rise immediately to a change 
of the spin background as a result of careful treatment of the constraint given in Sec. II B. 
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This is why the i-term is so efficient in destroying the AFLRO. Of course, there are many 
open questions in this approach, e.g., how to go beyond the MFA, what is the gauge field in 
this approach and what are its major effects and so on. These and other related issues are 
under investigation now. 



After submitting the original version of this paper, we found Ref. 56 where a similar 
approach was used to study the normal state properties of oxide superconductors. Apart 
from the difference in issues addressed in our paper and theirs a careful reader could easily 
discover the substantial distinction in the interpretation of transformations used in these 
two papers. 
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APPENDIX A: MATRIX REPRESENTATION OF THE PROJECTION 

OPERATOR 



The holon operators and h in the basis 



1 




1' 



•0' 

0^ 

1 



of holon states are given by 



while the spin raising and lowering operators S'^, S in the spin 1/2 space 
are given by 



1' 



1 




1 




In the product space \hole) ® \spin) the basis vectors are 

oj.^Voj„ 



ll,t) 







0" 



|o,T) 



0" 
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|o,i) 



'0' 



0" 
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/o\ 
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(Al) 



0" 



(A2) 



(A3) 



which form a complete set. 

The fermion-spin transformation defined by Eq. (7) gives the following matrix represen- 
tation for the fermion operators 

/O 0\ 



Ct = h^s- 



1' 




0^ 

1 



10 


Vo oy 
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0' 
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0' 

1 



/o 








1\ 


































0/ 


/o 








0\ 

















1 








\o 








0/ 


/o 








0\ 
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(A4) 



On the other hand, there are only three physical states in the constrained Hilbert space, 
namely 



(A5) 



To remove the extra degrees of freedom in the \hole) <^\spin) space, we introduce a projection 
operator P. By requiring P|l, t) = i) = |0), P|0, t) = | T), and P|0, i) = | i), one can 
easily find its matrix representation 

/I 1 0\ 





(l\ 




(0^ 






|0) = 





, IT) = 


1 


, U) = 


















and its hermitian conjugation 
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VO 1/ 



(A6) 



(AT) 



/I 0\ 
1 
1 

VO ly 

Using this projection operator, the electron operators in the restricted Hilbert space are 
given by 
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di ^Ph^s+p^ 



/O 1 0\ 



Vo oy 

/O 1\ 



vo 0/ 



/O 0\ 

1 

Vo oy 

/O 0\ 



VI oy 



(A8) 



as quoted in Eq. (11). It is then straightforward to verify the operator relations quoted in 
the main text Eqs. (12) and (13). In particular, the hole number operator 
/I 0\ 



rih 



1 



|0)(0| = -P(|iT)(iT| + |U)(U|)P^. 



(A9) 





vo oy 

The physical meaning of Eq. (A9) is transparent: The empty state should be counted 
only once, not twice. Since in the mean- field treatment the constraint on average doping 
concentration 5 is imposed directly on h^h, the sum rule for the physical electron is satisfied. 
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FIGURES 

FIG. 1. The momentum distribution of physical electrons in the mean-field approximation 

obtained by the fermion-spin transformation proposed in this paper (soUd line) in comparison 
with corresponding curves in the slave-boson (dashed line) and the slave-fermion (dot-dashed line) 
approaches. The doping concentration S = 0.5. 

FIG. 2. The spin background is assumed to be an antiferromagnetic state for the J > case 
of the ID t-J model. The t-J chain is squeezed and the spin configuration is rearranged due to 
the spin- up electron hopping: (a) Before hopping, when the holon is at site i. (b) After hopping, 
when the holon has hopped to site i + 1 from site i, while the spinon has hopped to site i from 
site i + 1. The spin polarization directions to the left of the hole are already optimized by the 
fermion-spin transformation, but there is still a hole in the chain, (c) Squeezing out the hole from 
the t-J chain. After this squeezing, the spin configuration is not optimal to favor the lowest spinon 
energy, (d) Rearranging the spin configuration from -oo to site i to favor the lowest spinon energy. 
The situation for the spin-down electron hopping is similar. 

FIG. 3. The momentum distribution of physical electrons in the mean-field approximation ob- 
tained by the fermion-spin transformation with the string fields (solid line), and without string fields 
(dot-dashed line), in comparison with the corresponding curve in the conventional CP^ approach 
(dashed line) (see Ref. 24). The doping concentration 6 = 0.25. 

FIG. 4. (a) The staggered magnetization, (b) the orbital current order parameter x, and (c) 
the order parameter ^ as a function of the hole concentration 6, for t/J=5 (solid line), t/J=3 
(dashed line). MF means magnetized flux phase with long-range order, while DF is the disordered 
flux phase. 

FIG. 5. Specific heat data as a function of temperature T (in units of J) at the hole concen- 
tration 6 = 0.2 for t/J = 3 (solid line), t/J = 5 (dashed line). 
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FIG. 6. The total energy Etotal (dashed hne) and the chemical potential /x/3 (solid line) as a 
function of doping S for t/J = 5. The range of the phase instability is roughly < S < 0.08. 
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